ON SOURCE-TYPE SOLUTIONS AND THE CAUCHY PROBLEM 
FOR A DOUBLY DEGENERATE SIXTH-ORDER THIN FILM 
EQUATION I. LOCAL OSCILLATORY PROPERTIES 



M. CHAVES AND V.A. GALAKTIONOV 

Abstract. As a key example, the sixth-order doubly degenerate parabolic equation 
from thin film theory, 

Ut = {\u\ m \u xxxxx \ n u xxxxx ) x in Ixl + , 

with two parameters, n > and m £ (— n, n + 2), is considered. In this first part of 
the research, various local properties of its particular travelling wave and source type 
solutions are studied. Most complete analytic results on oscillatory structures of these 
solutions of changing sign are obtained for m = 1 by an algebraic-geometric approach, 
with extension by continuity for m f» 1. 



1. Introduction: basic nonlinear model with complicated local and 

global properties of solutions 

1.1. Higher-order degenerate parabolic PDEs: no potential, monotone, order- 
preserving properties, not of divergence form, and no weak solutions. We con- 
sider the sixth-order parabolic equation from thin film theory, 

(1.1) ut = (\u\ m \u xxxxx \ n u xxxxx ) x in KxR + , 

with two parameters, n > and m G (—n, n + 2). This equation is doubly degenerate and 
contains the higher-order p-Laplacian nonlinearity |w aKE axna:| n (then p = 2 + n), and another 
one |ii| m of the porous medium type. The equation is written for solutions of changing sign, 
which is an intrinsic feature of the Cauchy problem (the CP) with bounded compactly 
supported initial data Uq{x) to be studied. 

The PDEs such as (II. ip . which are called sixth-order thin film equations (the TFEs- 
6), were introduced by King [IB] in 2001 among others for modelling of power- law fluids 
spreading on a horizontal substrate. Equation ( 11 .11) is quasilinear, where the diffusion-like 
operator includes two parameters m and n and is not potential (variational) and/or mono- 
tone in any functional setting and topology. It is not also an operator of fully divergence 
form, since the PDE admits just a single integration by parts, so that a standard defini- 
tion of weak solutions is entirely illusive. Of course, as a higher-order parabolic equation, 
(II .ip does not exhibit any order-preserving (via the Maximum Principle) features. 
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Overall, the TFE-6 (11 .11) . as a typical example of a variety of complicated nonlin- 
ear thin film models arising from modern applications, represents a serious challenge to 
general PDE theory of the twenty-first century, concerning principles and concepts of un- 
derstanding the common local and global features and properties of its solutions, which 
also need proper definitions. For other sixth-order TFEs including their derivation and 
mathematical properties, see a survey in [8], where further key references are traced out. 

It is well-known that nonnegative solutions of a wide class of higher-order TFEs can 
be obtained by special non-analytic and often "singular" e-regularizations and passing 
to the limit e — > + , that, in general, lead to free-boundary problems (FBPs). We refer 
to the pioneering work by Bernis and Friedman [3] and to the monograph on nonlinear 
parabolic PDEs [221 Ch. 4], where further references and results can be found. Actually, 
such singular e-regularizations, as e — > + , pose a kind of an "obstacle FBP", where 
the solutions are obliged to be non-negative by special free-boundary conditions that, in 
general, are not easy to detect rigorously. Solutions of the CP cannot be obtained by 
such techniques and require more involved and different analysis. In particular, analytic 
e-regularizations, with e — ► + , can be key for the CP, [8]. 

There is a large amount of pure, applied, and numerical mathematical literature devoted 
to existence, uniqueness, and various local and asymptotic properties of TFEs, especially, 
for the standard TFE-4: 

(1.2) u t = -(\u\ n u xxx ) x , where n > 0. 

Necessary key references on various results of modern TFE theory that are important 
for justifying principal regularity and other assumptions on solutions will be presented 
below and, in particular, are available in [T4"l[T5] and in a more recent paper [11]. See also 
[T3| Ch. 3], where further references are given and several evolution properties of TFEs 
(with absorption, included) are discussed. However, even for simpler pure TFEs such as 
fll.2p . questions of local and global properties of solutions of the CP, their oscillatory, and 
asymptotic behaviour are not completely well understood or proved in view of growing 
complexity of mathematics corresponding to higher-order degenerate parabolic flows. 

Using this, rather complicated, and even exotic, model equation (II. ip . we plan to ex- 
plain typical and unavoidable difficulties that appear even in the study of local properties 
of compactly supported solutions and their interfaces for higher-order degenerate nonlin- 
ear PDEs. We then intent to give insight and develop some general approaches, notions, 
and techniques, that are adequate and can be applied to a wide class of difficult degen- 
erate parabolic (and not only parabolic) PDEs. Since even the related ODEs of the fifth 
order for particular solutions get very complicated with a higher- dimensional phase space, 
we cannot rely on traditional ODE methods, which were very successful in the twenti- 
eth century for second-order ODEs, occurred for many particular self-similar and other 
solutions, with clear phase-planes. 

Instead, as a general idea, we propose to use parameter homotopy-continuity approaches 
using the fact that for some values of m, n (e.g., for m = 1 or m = n = 0, etc.), the ODEs 
can be solved by some algebraic-geometric methods, or leads to easier linear equations. 
Then, we use a stable "transversality-geometric" structure of the obtained solutions to 
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extend those into some surrounding parameter ranges. However, global extension of those 
solutions are not straightforward at all and often we are obliged to apply careful numerical 
methods to trace out some solution properties and their actual existence. 

Therefore, we are not restricted to quasilinear equations with semi-divergent operators. 
Without essential changes and hesitation, we may consider other fully nonlinear models 
such as the following formal parabolic PDE: 

(1.3) \ut\ a ut = (\u\ m \u xxxxx \ n u xxxxx ) x in lxR + (a > 0) , 

where a = leads to the quasilinear counterpart (11.11) . It is easy to propose other 
more artificial (and often awkward) versions of such PDEs without traces or remnants 
of monotone, potential, divergence form, etc., operators, which however can be used for 
applying our general mathematical concepts of analysis. 

1.2. On other models, results, and extensions. In the present first part, in Sections 
[2H7J we present a detailed study of some, mainly local oscillatory and sign- changing, 
properties of travelling wave (TW) and source-type solutions. In this connection, let us 
mention the first important pioneering results on oscillatory source-type solutions of the 
fourth-order quasilinear parabolic equation of porous medium type (the PME-4) in the 
fully divergence form with the monotone operator in H~ 2 : 

(1.4) u t = -{\u\ m - l u) xxxx in lxR + (m > 1), 

which were obtained by Bernis [1] and in Bernis-McLeod [3]; see also [221 § 4.2] for further 
details, and [9] for construction of a countable family of similarity solutions of (jl.4p . 

Notice that many key features of a local oscillatory structure of solutions of quasilinear 
degenerate PDEs with such operators do not essentially change not only for equations of 
higher, 2mth-order, but also for similar odd-order nonlinear dispersion equations (NDEs). 
In [S], as an illustration, we briefly review our approaches for the corresponding fifth-order 
counterpart of (ll.ip . which has the form (the NDE-5) 

(1-5) u t = {\u\ m \u xxxx \ n u xxxx ) x in Kx R + , 

Various odd-order PDEs occur in nonlinear dispersion theory. As a key feature, they ex- 
hibit finite interfaces, compacton behaviour, and shock/rarefaction waves; see first results 
in Rosenau-Hyman [19], a survey in [T3l Ch. 4], and [10] . as a more recent reference for 
shock wave behaviour for NDEs-5. The counterpart of (11.51) for m = n = is the linear 
dispersion PDE 

(1.6) u t = u xxxxx in Ix K + , 

which, by natural continuity-homotopy issues, gives some clues about oscillatory proper- 
ties of solutions for the nonlinear one (II. 5p . 

It is curious that the performed local analysis of oscillatory solutions close to finite 
interfaces remains similar for a number of other PDEs. For instance, in this sense, the 
parabolic equation (ll.ip has the counterpart which is a rather unusual PDE also belonging 
to the class of nonlinear dispersion equations 

(1.7) lift (1^1 \u xxxxx \ u xxxxx ^j 
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In its turn, a "local counterpart" of the nonlinear dispersion PDE (11.51) is the quasilinear 
sixth-order hyperbolic equation 



The second part of our paper [5] is devoted to a global construction of source-type 
solutions and some general aspects concerning the Cauchy problem for fll.ip . As a first 
general rule therein, we show that, for correct understanding the CP and some principal 
properties of its solutions, one should study the limit m,n — > + , where the linear in- 
harmonic equation occurs 

(1.9) Ut = u xxxxxx in Mx M + . 

One can easily construct for (11.91) the necessary TW and source-type (ZKB) solutions, 
to be compared with solutions of its nonlinear counterpart, the TFE ( II. ip . at least for 
sufficiently small m and n. As a principal feature, we observe that, for the CP, the 
solutions of both linear and nonlinear PDEs are oscillatory near finite interfaces (for 
( 11. 9ft . by implication, interfaces are at infinity, x = ±oo). 

In higher-order quasilinear degenerate evolution PDEs of parabolic, nonlinear dis- 
persion, or hyperbolic types, the identification of optimal local regularity properties of 
the FBP and the Cauchy problem becomes very difficult, to say nothing about general 
existence-uniqueness theory. Actually, we were surprised to observe that many classic 
techniques of PDE theory developed in the last fifty or so years not only cannot be ap- 
plied but also cannot provide us with necessary results in principle, since these are strictly 
oriented to lower-order PDEs. 

Evidently, related challenging difficulties appear even when constructing standard par- 
ticular self-similar solutions leading to higher-order ODEs. These cannot be studied by 
phase-plane analysis which is usual for the first or second-order equations. Moreover, 
several results, which have been obtained for given third or fourth-order PDEs by a com- 
plicated study of the topology of orbits, are often hardly applied to similar ODEs of the 
order that is higher by one. Many mathematical approaches are strictly attached to ODEs 
of the given and sufficiently low order. 

In this paper, we have specially chosen sixth- or fifth-order models ( 11.11) and (11.51) 
with complicated non-monotone and non-potential operators to describe a certain general 
scheme to study local and global properties of solutions including: 

(I) existence and uniqueness of travelling wave solutions, and 

(II) local oscillatory behaviour of solutions of changing sign near interfaces. 

In a forthcoming paper [5], we continue this study of (II. ip and will concentrate upon: 

(III) existence and uniqueness of source-type (ZKB-type) solutions, and 

(IV) general properties of solutions and proper setting of the Cauchy problem. 



(1.8) 




2. Travelling wave and source- type similarity solutions 



Let us now introduce two classes of particular solutions to be studied in detail. 
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2.1. Travelling waves. These are simplest solutions of nonlinear PDEs of the form 

(2.1) u TW {x,t) = f{y), y = x-\t (A G R), 
where, on substitution into (11.11) . / satisfies the autonomous ODE 

(2.2) A(/) = |/H/Wr/ (B) = -A/. 

This is obtained on integration once with a zero constant via the required flux continuity 
at the interfaces, where / = 0. 

For n = m = 0, the ODE (12. 2p is easy, with a linear bundle of exponential solutions: 

(2.3) = -Xf f(y) = e"», where /i 5 = -A. 

2.2. Source-type solutions. Equation (11.11) admits the standard self-similar solutions 

(2.4) u s (x,t)=r l3 F(y), y = x/t /3 , where (3 = ^^j, 
so (12. 4p preserves the total mass of the solution: 

Lus(x,t) dx = const, for all t > 0. 

Then F solves a slightly different ODE again obtained on integration once with the zero 
constant: 

(2.5) A(F) +{3Fy = |F| m |F (5) | n F (5) + (3Fy = in E. 

In view of the symmetry group of scalings, if F\(y) is a solution of (12. 5p . then 
(2-6) F a (y) = ^F 1 (^), with 7 = ^, 

is a solutions for any a > 0. Therefore, we are looking for a unit mass profile satisfying 

(2.7) jF(y)dy = l. 

As another normalization, one can take the condition -F(O) = 1. 

Remark on finite propagation. Finite propagation in the ODEs such as (12.21) and 
(12. 5p . written in a semilinear form (see below), i.e., existence of finite interface yo of 
arbitrary sufficiently small solutions, was well-known for a long time. We refer to the 
first ODE proofs in [U, H], [221 p. 392], and more involved energy estimates for general 
related higher-order elliptic and parabolic PDEs in [2l [20] and survey in [12] . However, 
we must admit that some extensions of energy methods to odd-order ODEs such as (12.21) 
or (12. 5p . i.e., for nonlinear dispersion operators involved, are not straightforward and can 
be technically rather involved. 

Thus, the principal question remained open is the oscillatory and non-oscillatory be- 
haviour of solutions close to interfaces. 
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2.3. Fundamental solution for n = m = 0. Then (12 .4p is the fundamental solution 

(2.8) b(x,t)=t-*F (y), y = x/t$, 

of the tri-harmonic equation (11.91) . The corresponding linear problem for the rescaled 
kernel denoted by F reads 

(2.9) F (5) + |F o2 / = in K, J F = 1. 

This has the unique solution F (y) by classic linear theory, j6]; see [5] for extra properties 
of the kernel F (y). Recall that we intend to use the linear rescaled kernel F Q given by 
(12. 9p in trying to understand the nonlinear one (12.51) . at least, for m, n ~ 0. 

3. Local existence of positive solutions and maximal regularity 

We now begin to study the behaviour of solutions of both ODEs (12 .2p and (12.51) near 
the interface point. For (I2.2p . we just assume that the interface is at y — and set 
f(y) = for y < 0. For (I2.5p . assuming that y < is the left-hand interface point of F, 
we perform the change y — y^ i— > y and obtain equation (12.21) with 

(3.1) A = (3yo < 0, 

up to an exponentially small perturbation as y — ► y . 

Thus, in both cases, we consider (12.21) in a neighbourhood of the interface at y — 0, 
where we scale out the parameter |A| 7^ 0, so that now A = ±1: 

(3.2) /W= T |/|«-i/ for y>0 } /(0) = (a = i=a 6 (-1, 1), A = ±l). 
Actually, the condition a G (—1, 1) holds in a wider parameter range 

(3.3) n > -1, mG(-n,n + 2). 

The assumption a > — 1 (m < n + 2) is purely technical that simplifies local analysis 
of f(y) close to "transversal" zeros and makes it quite standard. For a < —1, some 
technicalities occur that are often not in the focus of the present study. 

Recall that, in (13. 2p . both A = ±1 can occur for the TWs, and, always, A = — 1 for the 
source-type solutions. 

We should clarify what kind of solutions we are looking for, and, namely, which extra 
conditions are supposed to be posed at the interface y — to create proper solutions of 
the CP. For instance, one can look for solutions / 6 C* 4 ((— 6, 5)), with a constant 5 > 0, 
so this demands 

(3.4) /(0) = /'(0) = f (0) = AO) = / (4) (0) = 0, 

but it is not clear whether these correspond to the actual regularity associated with the 
CP. Bearing in mind more sophisticated ODEs that occur from models like (ll.3p . we would 
like to avoid any weak definitions of solutions that are based on integration by parts (for 
instance, for fully nonlinear operators, this makes no sense). 
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As a hint, we can compare the desired regularity for the CP with that for the standard 
zero height, zero-contact angle, zero curvature, and zero-flux FBP, for which the free- 
boundary conditions take the form 

(3.5) /(0) = /'(0) = /"(0) = and (\f\ m \f®\ n f®)(0) = 0. 

Therefore, at the interface y = 0, solutions of the FBP satisfy for any C > and some 
Ci = Ci(C7)eM, 

(3.6) /(y) = Cy 3 + W + - for y > / G C 2 ((-<5, 5)) (FBP), 
which is true provided that 

(3.7) 7 = 8+5n-3m > 3 ^ m < 5^n_ 

Of course, (13.61) exhibits less regularity at y = than (13.41) . In fact, the solutions of the 
FBP do not need and/or admit the zero extension for y < (by the definition, this makes 
no sense for the FBP). If they do, this already corresponds to the CP. 

For the CP, we need to use the concept of the maximal regularity, which for the ODE 
(13.2j) simply means that setting /(y) = for y < yields 

(3.8) solutions that are maximally smooth at y = admitted by the ODE (CP). 

The actual maximal regularity associated with the ODE under consideration needs special 
local analysis close to y = + . It is worth mentioning now that ( 13 Ah are not conditions 
of maximal regularity in general, for arbitrary values of m and n. 

3.1. Existence, uniqueness, and nonexistence of positive TW solutions. We be- 
gin by noting that, for A = —1, equation (I3.2p admits the positive solution 

(3-9) fo(y)=¥>oy fl , where ^=^ = ^1 and 

n + l 

( 3 - 10 ) ^= L (At _ 1)( ,_4_ 3 )(,-4)) ] m+ " (A = -l). 

This formula makes sense in the intervals /i G (4, +oo), (2, 3), and (0, 1), which give some 
relations between parameters m and n. For A = +1, two intervals fi G (1,2) and (3,4) 
are accepted. Note that the regularity of the solution (13.91) at y = (with the trivial 
extension f(y) = for y < 0) increases without bound as a — > 1~, since 

(3.11) /i = — > +oo as a — » 1~. 

For instance, it follows that fo G C 4 at the interface provided that 

(3.12) m<^ (/iG (4, +oo)). 

Below, we will explain the meaning of these positive solutions in the Cauchy problem. 

Notice that, for A = +1, in the most convenient interval (I3.12p . the ODE (13.21) does 
not admit solutions that are strictly positive in an arbitrarily small neighbourhood of the 
interface at y = 0. 

We first consider the natural "smooth" version of solutions of the maximal regularity, 
where the conditions (13. 4p hold, so that /o(y) is such a solution in the interval (I3.12p . 
The following classification of possible solutions holds. By oscillatory solutions of (13. 2p at 



the origin, we mean those that have infinitely many sign changes (for instance, isolated 
transversal zeros) in any arbitrarily small neighbourhood (0, 5) of the interface at y — 0. 



len: 



Proposition 3.1. Let ( QHD hold. The 

(i) for A = +1, any f(y) ^ for y « + of the problem (\3.2§ . ( |ff.^P is oscillatory at 
y = + , and 

(ii) for A = —1, there exists a unique local strictly monotone and positive solution of ( EH)) , 
which is given by ( \3.9f) . and all other solutions are oscillatory as y —>■ + . 

Remark: first comparison with the linear ODE. The existence of a ID manifold of 
positive solutions (with the interface parameter y ) is in good agreement with the same 
conclusion for the linear equation (12.31) . where the ID manifold of positive exponentially 
decaying functions at the "infinite interface" as y —>■ — oo is given by 

/(!/)= Ce», C>0. 

The rest of decaying orbits are oscillatory as y — > — oo. We thus prove that a similar 
property persists for some m, n > in the interval (I3.12p . 

Proof of Proposition l3~l[ (i) We argue by contradiction. Assume that f(y) is a nontrivial 
solution of (13.21) with A = +1. Then, if f(y) is positive in an interval (0,y*), it follows 
from the equation (13.21) that f^ < for < y < y*. After integration and taking into 
account the conditions (13. 4p . we obtain that f'(y) < and f(y) < on (0, y*), whence 
comes the contradiction. Similarly, there is no a negative solution on (0,?/*). 

(ii) The key idea of the proof relies on the translation invariance property of the equa- 
tion. For A = — 1 and positive solutions, consider equation (I3.2p in the form, 

(3.13) /(8) = |/|«-i/ for y>0 ( a = ^<l). 

As a first consequence, one can see after integrating this equation that any positive solution 
f(y) of the ODE problem (I3.2p is also convex and strictly increasing. 

We next prove that, for A = +1, the function f = fo(y) defined in (13 . 2 [) — (13.91) . is the 
unique positive solution of the problem (13.21) . (13. 4p . We begin by showing that if fo(y) 
and f(y) are different positive solutions of fl 3 . 2 j) . then either they are ordered, or they 
have mutually oscillatory behavior close to the origin (in the sense that they intersect 
each other infinitely many times at different values of y approaching the origin, i.e., the 
difference fo(y) = f(y) is oscillatory at y = 0). Notice that, in any domain of analyticity 
of f(y) > (where f'(y) ^ 0), all intersections with the analytic fo(y) are isolated points. 

We begin with the case a > 0. Assume that the assertion is false and that there exists 
y* > such that fo(y) and f\y) are different functions ordered in (0,y*) and such that 
fo(y*) = f(y*)- If fo(y) < f(y) in this interval, it readily follows from the equation in the 
new form (13. 13[) and the regularity properties of the solutions that fo 5 \y) < f^(y) for 
< y < y*. Hence, we obtain by integration in (0, y*) that fo(y*) < f(y*), a contradiction 
with the assumption. The same argument applies to the opposite inequality and therefore, 
the assertion holds. 



In order to establish the uniqueness result in both cases (ordered solutions or with mu- 
tually oscillatory behavior), we apply similar arguments to auxiliary solutions constructed 
by using the invariance translation properties of the equation. Assume for contradiction 
that there exists yo such that fo(y) ^ f(y) with /o(?/o) < f(yo)- Consider the auxiliary 
solution of (13. 2p . (13.41) defined as follows: 



where e > is chosen small enough such that fo(yo) < g(yo)- So defined, we have 
that g(y) is a non-negative solution of (13. 2p and, clearly, it satisfies by construction that 
g(y) < f(y) fory>e close enough to e. Therefore we have by continuity and the previous 
arguments that there exists a value y* G (e,y ) such that g(y) < f(y) in the interval 
{0,y*) and g{y*) = f(y*)- Hence, it follows from equation (I3.13P that ftf\y) < g^{y) for 
< y < y* and a contradiction at y* follows after integrating this inequality as above. A 
similar contradiction argument applies if we assume that /o(z/o) > f(yo) by considering 
the auxiliary function g{y) = f (y — e). 

The uniqueness result for a < is obtained by using similar arguments but taking into 
account that in this case h(s) = is a decreasing function for s > 0. This fact gives 

after subtracting and integrating five times that 



whence if the solutions are ordered close to the origin, the signs of the left- and right- 
hand sides of this expression are different, so the contradiction readily follows. The same 
argument applies to g(y) constructed as above, so mutually oscillatory behavior is also 
disregarded in this case. This completes the proof. □ 

3.2. Positive solutions by fixed point theorem. Indeed, we were lucky to have the 
explicit positive solution (13.91) . which is connected with the invariant scaling group of 
transformations of the ODE (13.21) . We now sketch another approach to detecting a unique 
positive solution for more arbitrary nonlinearities without using and relying on explicit 
calculus via a scaling group. 

We again consider the ODE (13.131) with A = —1, i.e., look for positive solutions of 
(13.131) . where f a can be replaced by more general functions q(f) > for / > 0. The main 
non-uniqueness difficulty for (13.131) is that it always admits the trivial solution f(y) = 
if a > 0. However, this disappears for the inverse function y = y(f), for which 




for y > e, 
for y < e, 



y s 



M - fo(y) = J - J(f a - fo a )(r) dr ... dz, 







£i J_ ell f_l_V— f 

J — yll J — \ yl j yll J 
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Then we obtain the following ODE for y(f): 



(3.14) 



— (((((?/) y>) y>) y>) y') 



f a for / > 0. 
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Multiplying by y' and integrating over (0, /) with the zero boundary condition yields 

((((7)7)7)7)7 = //v,d/ = ry -ajr-tydf. 

Multiplying again by y' and integrating gives on the right-hand side a smooth quadratic 
integral operator depending on y 2 , etc. 

After five integrations like that, we obtain an integral equation for y(f) of the form: 

(3.15) y (f) = M(y)(f) for / > 0, 

with a smooth operator A4 of y containing the maximum fifth degree y 5 polynomial 
dependence on y. Therefore, this operator is a contraction in a space of continuous 
functions provided that the typical integrals such as j Q f a ~ l {-) df converge, i.e., 

(3.16) a>0 me (0,1). 

This restriction can be weakened by choosing special classes of functions y(f) with pre- 
scribed envelopes as y — ► 0. 

Finally, Banach's Contraction Principle (see e.g. [TT1 p. 206]) guarantees existence 
and uniqueness of a positive solution of f 1 3 . X 5 [) . and hence of (13.131) . For non-monotone 
changing sign solutions of (13.131) . this inverse function approach obviously fails. We will 
study such a behaviour in Section 15.11 by an extra scaling. 

3.3. Maximal regularity. It is easy to see that (I3.9P describes the best (maximal) 
regularity at the interface that is provided by the ODE (12.21) . Indeed, (13.91) established 
the only possible balance between two terms in (12. 2p . In other words, loosely speaking, if 
such a balance is violated, the solution must behave along the kernel (the null-manifold) 
of the nonlinear or linear operator on both sides of (12.21) . which is obviously trivial. Recall 
that this maximal regularity reflects the intrinsic properties of the CP. For the FBP, the 
regularity is different; cf. (13. 6p . 

Proposition 3.2. Let a > and let the solution o f (jff.gp . Qff.^D be written as 

(3.17) f(y)=y^(y) for y>0 ( fl = ^ = ^). 

Then <p(y) is uniformly bounded for y > 0. 

Proof. Let us prove that cp(y) is globally bounded, i.e., there exists a positive constant C 
such that 

(3.18) \<p(y)\<C for all y > 0. 
Assume that / ^ 0, and, for every y > 0, define 

(3.19) V = snp m {z G [0, y] : < \f(z)\ Vs G (0,2/)}. 

Using the ODE such as (13.131) and integrating five times over (0,y), we obtain by the 
definition in (I3.19P that 

\m\<c\f(y)\ a v 5 , 

so the desired inequality (I3.18P readily follows for y and, by (I3.19p . for every y > 0. □ 
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3.4. Nonexistence of positive source-type similarity profile. We finish this analysis 
by dealing with the question of nonexistence of positive source-type solutions. We prove 
that, in the basic interval ( 13. 12ft . the ODE problem ( 12. 5p has no positive solutions with 
the symmetry conditions at the origin. 

Proposition 3.3. There are no symmetric positive solutions with compact support of the 
problem (\2.5§ , with conditions (|ff.^P for F(y) at the interfaces. 



Proof. Let F(y) be a positive solution of (12.51) with compact support [—yo,yo\. It is clear 
from the equation ( 12.51) and the positivity assumption that F^ 5 >(y) > in (—y ,0). By 
integrating this inequality over (0, y ) and assuming the required regularity (13 .4p of F(y) 
at y — yo, we infer that F'(0) > 0. Hence the symmetry condition fails that completes 
the proof. □ 

Thus, the sufficiently smooth source-type profile in the CP must be oscillatory at the 
interfaces. We now begin to study the character of such oscillations. 

4. Existence and uniqueness for the initial value problem 

In order to analyze the existence and uniqueness problem, we introduce the following 
(cf. □): 

Definition 4.1. A function f is said to be a weak C A -solution of the equation (\2.2\i on 
an interval I, if f 6 C 4 and for any y , y e / , 

(4-1) f^{y) - fW{yo) = ±J \f(s)\«-if(s) ds. 

yo 

One can see that, for a > 0, solutions are C 5 and hence classical. For a < 0, if we use 
the asymptotics (13.171) with uniformly bounded (say periodic) oscillatory components (p, 
weak C 4 -solutions exist for 

(4.2) a>-\ (i.e., &>-l). 

Consider the Cauchy problem for the equation ( 12.21) with the initial conditions 

(4.3) fU)( yo )= aj , j = 0,1,2,3,4. 

It is clear that this initial value problem can also be written, by introducing the function 
h(y) = \f(y)\ a ~ 1 f(y)y i n the equivalent form, 

f^(y) = a,±h(y), 

h'(y) = \f(y)\ a - 1 f(y), 

/<%,)=«,•, J = 0,1, 2, 3, 
h{y ) = 0. 

From the standard theory of ordinary differential equations, for every a > 0, problem 
(I4.4p is solvable in some neighborhood of y = yo. Let (/, h) be a solution. On the one 
hand, it is clear that, for a > 1, such solution is unique since the Lipschitz condition on 
the nonlinearity is satisfied. On the other, for a < 1, the only continuity is guaranteed 

n 



(4.4) 



and in fact uniqueness does not hold, for instance, if a>j = for every j = 0, ..,4, as it 
is shown by means of the construction of positive solutions in the previous section. A 
partial answer to uniqueness is established in the following: 



Proposition 4.2. Let < a < 1. If in (EO) 



EM>o, 

3=0 



then the initial value problem {4-4) has at most one solution in a neighbourhood ofy = ijq. 



Proof. The proof follows the ideas in [I], H] . We prove uniqueness in a small right-hand 
neighbourhood of y = yo- Uniqueness on the left is obtained in a similar way. Denote by 
j* the smallest value of j such that aj is non-trivial, and, for such j = j*, introduce the 
function 

/ \ f / fiiv \i , if y > yo, 
l9i{yo) = 0, if y = yo- 
Assume for contradiction that there exist two different solutions fi and f 2 for small 
y — yo > 0. Due to the invariance properties of the equation, we may assume without 
loose of generality that fi > f 2 > 0. It is clear by integrating the equation and taking 
into account the initial conditions that 

\A(y) - h{y)\ <J...J\\fi\ a - l fi(r) - \f 2 \ a - 1 f 2 (r)\dr...dz. 
yo yo 

Hence, after dividing by (y — yoY , taking into account that h(s) is Lipschitz continuous 
away from s = and integrating five times, we get 

\9i{y)-M\ <C{y-y Q ) A+{1 - a)j max \ 9l (s) - g 2 (s)\. 

[yo,y] 

Since 4 + (1 — a)j > for any j < 4, the contradiction follows for y — y > sufficiently 
small. The proof is complete. □ 

Next we introduce a comparison principle which is valid for a > with A = — 1 (the 
positive sign on the right-hand side in (13.21) ) and a < with A = 1. 

Proposition 4.3. Let fx and f 2 be two solutions of equation fl^.ip on [y , oo) satisfying 



fi\yo) ^ /s \ (yo) f or 3 — 0, 1, 2, 3, 4, with strict inequality for at least one of them. Then, 
h{y) > f2(y) for every y > y . 

Proof. The inequality is obvious by continuity and the assumptions on the initial condi- 
tions on an small interval [yo,yo + £)■ We see that in fact the inequality holds for every 
y > yo. If the assertion is false, we denote by y* = sup{y < yo : fi(s) > f 2 (s), s G 
(yo,y)}- It is clear that, at this point, fi(y*) = f2(y*)- However, for a > and A = —1, 
we obtain after subtracting and integrating five times, 

h{y*) - / 2 (y*) = E (a? - a{) + J ... JdM^Mr) - |/ 2 |— Va(r)) dr ... dy > 0, 
3=0 yo yo 
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leading to a contradiction. In a similar way, a contradiction is obtained for a < and 
A = l. □ 



5. Oscillatory component: dissipative systems, periodic behaviour, and 

heteroclinic bifurcations 

5.1. Dissipative dynamical system for the oscillatory component. We now study 
solutions of (13. 2p of changing sign. Bearing in mind (13.171) . we will be looking for solutions 
in the form 



(5.1) f(y) = y"<f{s), s = \ny, where \i 



5 _ 5(n+l) 
I— a m+n 



Here if is called the oscillatory component of the solution, and is a sufficiently smooth 
function; see below. The positive solution (13.91) corresponds to the particular case, where 
ip(s) = ip (see (13.101) ). and this exists in the corresponding parameter ranges. Besides 
( 13 .91) . there exist many other changing sign solutions. 

According to Proposition 13.21 we know that, at least for a > 0, the oscillatory compo- 
nent satisfies 

(5.2) (p(s) is bounded and y4 as s — ► — oo. 

Therefore, (15.11) gives a clear picture of typical regularity of solutions at the interfaces, 



/ G C 4 ((—5, 5)) for m< 



n+5 
4 ' 



(5.3) / G C 3 ((-5, 5)) for =±5 < m < ^±5, 

/ G C 2 {{-5, 5)) for 2r*±5 < m < 3^f5. 

According to (13. 7p . in the first two ranges the CP always has a better regularity than the 
FBP. 

Substituting the representation (15.11) into the ODE (13.21) . after simple manipulations 
via scaling properties, yields the following autonomous equation for (p(s): 

(5.4) P 5 (ip) = TM°~V m R (A = ±l), 

where a = jt^- Here P5 is an easily derived linear differential operator of the form 

p 5 ((p) = e(-^- 1 ) s (e- s (e- s (e- s (e- s (e /i V) , ) , ) / ) / )' 
= y^ (5) + 5(/i - 2)tpW + 5(2^ 2 - 8n + 7)ip"' 
^ ' ' + 5(/x - 2)(2/i 2 - 8/x + 5 V" + (5/i 4 - 40/i 3 + 105/i 2 

- 100// + 24V + - l)(/x - 2)(/i - 3)(/i - 4)(^. 

The best and simplest connection satisfying (15. 2p with the interface at s = —00 is 

(5.6) <f = v?*(s) is a (non-constant) periodic solution of (15. 4p . 

On the other hand, any uniformly bounded solution of (15. 4p for y <^ — 1 will fit. 

We now concentrate on periodic connections. We list the following properties that 
approach us to existence of a periodic orbit of changing sign: 
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Proposition 5.1. Let a G [0, 1) and (cf. ([Hp ) 

(5-7) ^=I^ = ^>4, i-e., m<n±5. 

T/ien the fifth-order dynamical system ( 5.4\) satisfies: 

(i) no orbits are attracted to infinity as s — > +oo ; and 

(ii) z£ is a dissipative system with bounded absorbing sets defined for any 5 > by 

(5.8) limsup s _ +00 \(p(s)\ G B* s = [-C* - 5, C» + 5], w/iere C* = (5!)~^. 

The absorbing set in M 5 includes all vectors (tp, tp' , tp" , tp'" , Lp(^) T , where all derivatives 
are uniformly bounded by a certain constant. 

Proof, (i) The operator in (15.41) is asymptotically linear [TTl p. 77] with the derivative at 
the point at infinity P§{(p) that has the characteristic equation obtained by setting 

(5.9) V = e ps =^ p 5 (p) = (p + p)(p + p-l)(p + p-2)(p + p-3)(p + p-A) = 0. 
Therefore, all eigenvalues are real negative or non-positive: 

(5.10) p k = k-p<0 for k = 0, 1,2,3,4, provided that 4 - p < 0, or m < 2±5. 

Thus, <p = oo cannot attract orbits. 

(ii) Actually, this follows from (i) by using an extra scaling. On the other hand, this is 
easy to see from the original equation such as (13 . 1 3[) for A = —1 (for A = +1 the proof is 
precisely the same). Taking an arbitrary solution f(y) defined for all y > 0, we integrate 
(13.131) five times to get that 

f(y) = sup z6(0ia) \f{z)\ > 

satisfies the inequality 

/K11 , f(y) = C + Cry + CV + C 3 y 3 + C 4 y 4 + J J J J J \f(z) \«- l f{z) dz 

(5.11) o 

< Co + C x y + C 2 y 2 + C 3 y 3 + + i y 5 f a (y), 
where Ci > are some constants. It follows that 

(5.12) f(y)<C + C.yT^ for all y > 0, 
provided that > 4. In the variable (15. ip . this yields (I5.8p . □ 

5.2. Periodic oscillatory component. In connection with Proposition 15.11 existence 
of a periodic orbit for dissipative dynamical systems is a standard result of degree theory; 
see [T71 p. 235]. Nevertheless, classic theory deals with non-autonomous periodic systems, 
with the given fixed period. Therefore, these results do not directly apply to the fifth- 
order dissipative autonomous dynamical system (15.41) . so we need some extra arguments 
to justify existence of periodic solutions for some values of a and how this disappears 
when a gets negative. 

It turns out that it is easier to establish first where such periodic orbits are nonexistent. 
The existence evidence will be presented later on. 
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Coefficients of the polynomial operator PM) 




FIGURE 1. Coefficients of the polynomial operator P§{<p) in (15, 5j) and (15.151) . 

Proposition 5.2. (i) The ODE ( 5.4 ) with A = — 1 does not admit a nontrivial T^-periodic 
solution (p* for 

(5.13) ne(^i4) or a e (-1,-0.9655...), where ^ « 1.45608... , ^ « 2.5439... . 

(ii) Analogously, the ODE {5.4) with A = +1 does not admit a nontrivial T^-periodic 
solution <p # for 

(5.14) 5 



3.22474.. 



2 < fJ. < fJ. 
Proof. We write down (15.41) as follows: 
(5.15) P 5 (<p) 



or aE (-l,a+), 



where oil 



fj,* — 5 



— 0.5505. 



V? (5) + a 4 <y9 (4) + a 3 ip'" + a 2 ip" + a\tp' + a ip 



Ti^r-v (a=±i), 



where the coefficients = aj(/i) are as in (15.51) . For convenience, we present graphs of 
these coefficients in Figure [H 

Multiplying (I5.15P by ip and <p' in L 2 (0,T*) and integrating by parts yields, 

(5-16) a 4 J(if'l) 2 - a 2 jVJ 2 + a J = T J \<p m \°+\ 

(5-17) /M 2 -as/(^) 2 + ai/(^) 2 = 0, 

respectively. It is not difficult to see that those identities make sense for all a > —1, so 
this restriction enters all our conditions. 

Note that multiplication by <p" gives a new formal identity 

(5.i8) -a* jv;t + « 2 jko 2 - « jvj 2 = ±a / i</>*r v*) 2 - 

However, on the right-hand side, there appears, after integration by parts, an integral 
with unknown converging properties, so using it could lead to wrong conclusions. More 
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precisely, this integral makes sense for a > only, which stays far away from the parameter 
ranges of interest. 

Thus, here and later on various integrals of ip* and its derivatives appear. The con- 
vergence of these integrals for a < is checked by using local properties of solutions of 
(15. 4p . In particular, these are easy for a > — 1, which justifies identities (15.161) . (I5.17P in 
the range of parameters under consideration, (I3.2p - (I3.3I) . Such convergence can be quite 
tricky for negative a < — 1 that are not under scrutiny in the present study. 

(i) To begin with our nonexistence purposes, it suffices to use just the single second 
identity (15.171) . This gives that ip* ^ const, is nonexistent if 

(5.19) ai > and a 3 < 0. 

It is easy to prove that there exist values and [if roots of a± such that a\ > and 
a 3 < on and hence, (p* does not exist, at least, on such interval. Numerical 

approximation of these roots shows that « 1.45608... and fi^ ~ 2.5439... . The results 
are illustrated in Figure [H 

(ii) In a similar way, the nonexistence result for A = 1 is extended to (I5.14p . by using 
identity (15.161) . This directly gives that tp* ^ const, does not exist if 

(5.20) a 4 > and a 2 < 0, 

which is true if /i G (§,/^!~], where /i+ = 3.22474.. is the largest root of a 2 . □ 

We continue with establishing the principal fact that the periodic solution ip* is hy- 
perbolic in some parameter ranges of interest. This establishes an important corollary 
concerning the type of heteroclinic bifurcation, and where it appears from; see below. 
This also reflects the stability and instability properties of periodic solutions for A = +1 
and A = — 1. We concentrate on the more important case A = — 1, which in (15. 4p corre- 
sponds to oscillatory behaviour of source-type solutions. 

Proposition 5.3. Let <p*{s) be a non- constant T* -periodic solution of (5.4), A = —1, for 

(5.21) /iG(/4,3], i.e., a G (-0.9655.. .,-§]. 
Then tp* is hyperbolic. 

Proof. For convenience, we write down (15. 4p . A = — 1 as 

(5.22) <p*: A^( V ) = -P 5 ( V ) + \^\ a - 1 V = 0. 
Consider the eigenvalue problem for the corresponding linearized operator 

(5.23) A.L(<P*)1> = -AW - M = A fc ^, 

where we use that a < 0. Note that the potential here can be rather singular at zeros 
of y?*(2/) that can essentially affect the setting of the eigenvalue problem and hence the 
actual regularity of eigenfunctions. In particular, we always assume that these zeros are 
transversal, so that the potential |<^*(2/) I"" 1 at zero, say, at y = is not singular as ~ 
as y — > 0. Then by Hardy-Rellich-type inequalities, we have got necessary embeddings 
that guarantee compactness of the resolvent and hence discreteness of the spectrum, which 
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allows to deal with the hyperbolicity issue for the periodic solutions. In further simple 
manipulations, we naturally assume that all these are justified. 

Assume first, for simplicity, that an eigenvalue A& is real. Multiplying ( 15. 23ft by ip and 
integrating over the period (0,T*) by parts yields 

(5.24) -04 jm 2 + «2 jm 2 -a j^ 2 -\a\j \^\ a ~^ 2 = X k j ^ '. 
It follows that 

(5.25) Afe < (ImAfc = 0), provided that a 4 > 0, a 2 < 0, a > 0. 

This gives the //-interval in ( 15.211) . For A& G C, similarly, first multiplying in L 2 by the 
complex conjugate ip, and next by ip the complex conjugate ODE, summing up yields 

(5.26) -04 / |VT + faj \f\ 2 -a J |^| 2 - \a\ j I^MVf = ^J \M 2 < 0, 
i.e., ReAfc < 0, under the same inequalities for the parameters as in (15.251) . □ 

The identities ( 15.241) and (15.261) can be used to detect other hyperbolicity (and stability) 
parameter ranges of and also for A = +1. 

As an important corollary, we obtain that, by classic bifurcation-branching theory and 
implicit function theorem for periodic solutions [2TJ Ch. 6], for A = —1, a hyperbolic 
periodic solution <p*(s), existing at some parameter value /i G 3), can be extended 
into a small open neighbourhood of (/io — fa + This implies the following: 

Corollary 5.4. The ^-parameter domain of existence of aperiodic solution tp*{s) of '(5.4), 
A = — 1 for /j, G (2,3), if it is not empty, contains a connected interval (/Ah, fa), where 
fa > fJ>t and fa < 3. 

Finally, we expect that the periodic connection (15. 6p together with its ID stable mani- 
fold as s — > — oo is the only transition to the interface point at s = — oo (y = + ), though 
this is difficult to justify completely rigorously. Also, we cannot prove uniqueness (see 
below) of a periodic solution of (15. 4p . Before, the only result on existence of a periodic 
orbit is obtained in (7J § 7.2] for the third-order ODE like (15.41) with n = and replacing 
the linear operator, 

P 5 (cp) ^ P 3 (<p) = if'" + 3(> - iy + (3/i 2 - 6/i + 2)<p? + - l)(/x - 2)<p, 

by using the fact that this dynamical system is dissipative. Then classic theory [TTl § 39.3] 
applies together with a shooting- type argument. This results was improved in similar lines 
in [TT| § 5] with a sharper estimate on the periodic solution existence interval: 

(5.27) < n < n h G (§,n+), where n + = = 1.9019238... . 
The actual heteroclinic bifurcation occurs at 

n h = 1.7599... (the TFE-4), 



and was calculated numerically, [7]. 
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5.3. Heteroclinic bifurcations of periodic solutions. The ODE (j5.4p is fifth-order 
quasilinear, is rather non-standard, and has a non-Lipschitz nonlinearity. Moreover, pe- 
riodic solutions exist not for all m,n > 0, namely, on some open (cf. Corollary 15.41) 
interval m G (0, m^n)), with the bifurcation exponent to be discussed. The classic 
results on existence of nontrivial periodic solutions based on rotation vector field theory 
[T7J p. 50] or branching theory [2TJ Ch. 6] do not apply to such ODEs. Solutions that are 
bounded on the whole axis (which is fine for (15. 2p ) also cannot be detected along the lines 
of classic theory in [T71 p. 56]. Moreover, we are interested in solutions that are bounded 
(with non-zero limits) as s — > — oo only. Several other techniques for existence of periodic 
solutions also fail; see references and comments in [T3J p. 140]. Therefore, we will rely 
on careful numerical evidence, especially when talking about the uniqueness of periodic 
orbits and their stability. 

The TW case A = +1 . In Figure [2} we show a stable periodic behaviour for the ODE 
(15 .4p for m = 1 and n = (so this is the standard thin film case) for A = +1. It turns out 
that such a periodic solution persists until a heteroclinic-like bifurcation which occurs at 

(5.28) m h = 1.337968147... (n = 0, A = +1). 

To compare with Proposition 15.21 we present the corresponding universal critical values 

(5.29) = = 3.7370... G (3,4), a h = ^ = -0.3380... . 

In terms of the original parameters m and arbitrary n, we have that the heteroclinic 
bifurcation occurs at 

(5.30) m h {n) = JL + (JL - l)n « 3.7370 + 2.3737n. 

Note that, for A = +1, (15.301) is essentially far from the predicted nonexistence value 
(15.131) . It is interesting to check whether the system (15.161) . (I5.17P (plus other identities 
if any) can supply us in this case with a better estimate of /ih- 

This scenario of a heteroclinic bifurcation of stable and hyperbolic periodic solutions 
is typical for dynamical systems; see [HI Ch. 4]. Indeed, the hyperbolicity property as 
in Proposition 15.31 excludes saddle-node bifurcations of periodic solutions (these demand 
existence of a A G i R) , at which <p* can disappear. Therefore, the appearance of a (stable) 
heteroclinic orbit as /i — > /ih in such dynamical systems is most plausible. 

Figure [3] shows this typical formation of a heteroclinic orbit — y?o — > V^o as m ~ > m h- 
Here ±y?o are constant equilibria of (15 ,4p given as in (I3.10p by 

n + l 

( 5 - 31 ) ^°-[- M,-D(.-4-3)(.-4) ] m+ " ( A = +!)- 

Since by ( 15.29j) . for n = and m = m^, we have ^ G (3,4), both equilibria ±y?o exist. 
These results have been obtained by the MatLab (the ode45 solver) with the enhanced 
accuracy parameters Tols= 10~ n and the same parameter of regularization in both de- 
generate and singular terms in (15.41) . 

Notice that, in view of (13.101) . solutions (15.11) for m G (|,mh) do not satisfy the last 
condition in (13.41) . i.e., 

f^\y) is not bounded as y — > + . 
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FIGURE 2. Stable periodic behaviour for the ODE (pT3|) for m = 1, n = 0, A = +1. 

Nevertheless, the structure (15 .11) corresponds to the maximal regularity for the ODE (13.21) 
and correctly establishes the balance of its nonlinear and linear terms. In Figure HI we 
show a stable periodic behaviour for m = n = 1 and A = +1. 

The source- type case A = -1 . The case A = -1 that includes the source-type solu- 

In this case, the periodic behaviour is unstable 



tions, is shown in Figure \5\ for m — n — 1 
as s — > +oo (as well as s — > — oo), though is clearly visible. For n 
bifurcation occurs at (see details in [131 P- 142]) 

(5.32) m h = 1.909... (n = 0, A = -1). 

The corresponding other critical values are 



0, a "heteroclinic" 



(5.33) 



2.619... e (2,3) and a h 



-0.909. 



l-«h ' Mh 

Observe that this «h is sufficiently close to the nonexistence one a* = —0.9655... in ( 15.131) 
in Proposition 15.21 The corresponding critical values of mh(n) are then given by 

(5.34) m h (n) = ^ + - l)n « 1.909 + 0.909 n. 

Let us state the following conjecture on periodic solutions and their stable manifolds. 
Conjecture [5l2. For all n > 0, there exists m^(n) given by < \5.30b for A = +1 (resp. by 
(5.34) f or ^ = — 1); 5mc/i that for all m e (0,mh(n)), t/ie OL'ii' (5.^) wi/i A = +1 (resp. 
A = -l); 

(i) /ias a unique stable (unstable) as s — ► +oo periodic solution <^*(s) of changing sign; 

(ii) as s — > — oo, t/ie periodic solution (p*(s) is unstable and has a ID stable manifold; and 

(iii) t/iere exists a heteroclinic bifurcation atm^(n), so ip*(s) is nonexistent for m > mh(n). 
It is worth mentioning the obvious instability of (p*(s) as s — > — oo (while approaching 

the interface point at some y = yo) is associated with the translational invariance of the 
ODEs such as (I3.2p admitting shifting in yo. The precise meaning and the significance of 
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FIGURE 3. Formation of a heteroclinic orbit asm-* = 1.3380... for n = 0, 
A = +l. 

m=n=1, X=+1 
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FIGURE 4. Stable periodic behaviour of (|5.4p for m = n = 1, A = +1. 

the conclusion (ii) will be explained in the second part of the paper [5] , which will be key 
for a well-posed shooting of a global source-type similarity profile. 

5.4. Oscillatory component: nonlinear dispersion model. Substituting TW solu- 
tions (12. ip into fll.7p and integrating twice with zero constants of integration yields the 
same ODE (12. 2p with the only change on the right-hand side: 

-A i-> A 2 > 0. 

Therefore, for the oscillatory component we obtain the same equation (15. 4p . and in ex- 
amples and figures presented there, we always take A = — 1, so that —A = 1 > 0. In 
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FIGURE 5. The trace of an unstable periodic orbit of (|5.4p for m = n = 1, A = —1. 

particular, the oscillatory behaviour near interfaces persists until the bifurcation expo- 
nent (15321) for n = 0. 

6. Existence and uniqueness of oscillatory TW solutions for m — 1 

In this section, we analyze the problem of existence and uniqueness of the oscillatory 
travelling waves solutions, as well as the character of their oscillations. We begin by 
proving existence of the oscillatory TW solutions via a different approach. 

6.1. Existence and uniqueness: an algebraic-geometric approach to periodic 
orbits. Here, we develop an alternative approach to existence of periodic orbits of (15.41) 
for both A, positive and negative. We focus our attention on the construction of solutions 
with right-hand side interfaces, f + (y) and f"(y) corresponding to A > and A < 
respectively. By reflection, we obtain solutions with left-hand side interfaces, namely 
f + (—y) for A < and f~{—y) for A > 0. 

We begin with the original equation for f(y), (I3.13p . which, for the case m = 1 (and 
putting for convenience |A| to make calculations easier), takes especially simple "weakly 
nonlinear" form 

(6.1) / (5) = 5! sign/. 

Let us explain the main ingredients of the strategy of our construction. 

By using the properties of the solutions, we may consider for convenience, that the 
positivity domain of f is given by (—1, 0), so /(0) = /(— 1) = and f(y) > on (—1, 0). 
Hence, such a function denoted now by fo(y) satisfies on (—1,0), 

(6.2) /( 5 ) = 5! f (y)=y(y + l)[a + by + cy 2 + y 3 }, 

where the constants a, b, c are chosen so that fo(y) > on (—1, 0). 
We next extend fo(y) to y > as follows in three steps: 
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(i) take -f (y), 

(ii) shift it to the right by y\ = 1 to get —fo(y — 1), and 

(iii) rescale it by the invariant scaling group for (16 .11) to get for y > 0, 

(6.3) fo(y) = — C 5 /o(g — l), with a parameter G > 0. 

In order to have a smooth solutions on [—5,5], we need four matching conditions on 
zero jumps of four consecutive derivatives, 

(6-4) [/a (0) = [/*] (0) = [ft"] (0) = [/ (4) ] (0) = 0. 

Thus, we obtain four algebraic equations with four parameters a, b, c, and G. They can 
be written as: 



(6.5) 



(1 - G 4 )a + G 4 b - G 4 c = -G 4 , 
(1 + G 3 )a + (1 - 2G 3 )b + 3G 3 c = AG 3 , 
(l + G 2 )6+(l-3G 2 )c=-6G 2 , 
(l + G)c = 4G- 1. 



We remark that by the argument of continuation of /o to a solution of the problem 
with a finite right-hand interface, the condition 

< G < 1 

is required. Therefore, we restrict ourselves to analyzing system (16. 5p in this range of 
parameters. 

We begin by noting that, for every fixed value of G G (0, 1), (16.51) is a linear system in 
a, b, c of four equations. Considering the three last equations, one can easily check that 
the determinant of the associated matrix 

d= (1 + G 3 )(1 + G 2 )(1 + G) 

is not zero, and, hence, for every fixed G G (0, 1), the system of the last three equations 
has a unique solution. Therefore, it will be the solution of the complete system (16.51) if 
the determinant of the matrix, 



(6.6) 



/l - G A G A -G 4 -G A \ 

l + G 3 1-2G 3 3G 3 AG 3 

l + G 2 1-3G 2 -6G 2 

\ 1 + G4G-1/ 



is zero. Denoting by D(G) the determinant of this matrix, we have that D(0) and D(l) 
are both negative. Hence, it readily follows that either no roots of D(G) are available, or 
at least two roots exist. After analyzing in more detail D(G), one can see that there exist 
two roots, that can be obtained numerically: 



Gi = 0.178318... and G 2 = 0.7060378... . 
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However, the positivity condition of / in (—1,0) implies a < or, equivalently, by using 
the first equation in (I6.5p . gives 1 + b — c > 0. Manipulating the two last equations in 
( 16.51) . one obtains that this restriction holds, if only if, 

F(G) = 3G 2 - 10G + 3 > 0. 

One can check that this is satisfied by G\ and allows to disregard G 2 . 

It is key that the existence of the value G 2 also plays an important role in the analysis, 
since it provides, by means of the construction explained above, solutions of the equation 
with the left-hand interface for A < 0, i.e., 

(6.7) / {5) = -5! sign/. 

These arguments allow to state the following: 

Theorem 6.1. (i) The above algebraic system (I 6. <5|) has a unique solution Gi E (0, 1) 
satisfying the positivity condition of f = fo(y;Gi) in (—1,0) and a unique solution G 2 E 
(0,1) satisfying the negativity condition of fo = fo(y;G 2 ) in (—1,0). 

(ii) Equation ( 5.$ , with the sign " + " {respectively " — ") and a = 0, has a nontrivial 



periodic solution f + (y) (respectively f (y))> with the period 
(6.8) T* = 2 In G\ {resp. T* = 2 In G 2 ) ■ 

Proof, (i) We prove the result for A positive and the solution corresponding to G\. For 
A < the same arguments apply. 

Once a solution of the algebraic system (16.51) has been found, we can extend the solution 
fo(y) similarly to the next interval of oscillations, etc. indefinitely. Then the sequence 
of positive and negative humps will converge according to some geometric series. For 
instance, the right-hand interface of the interval of every extension, given by the standard 
geometric series, 

b n = G + G 2 + G 3 + ... + G n (G = G x or G 2 ), 

converges as n —>• oo to the right-hand interface of the solution, ^-y. 

It is not hard to see that, after reflecting and moving the interface to y = + , this gives 
precisely the behavior (15.11) . i.e., 

(6-9) f (y) = (y -y) 5 ^(s), s = \n(y -y), 

where j/o = denotes the right-hand interface of fo. We next prove that the oscillatory 
component </?* is periodic with the period (16. 8ft . In fact, we show that satisfies a 
stronger property, 

Lp*(s + lnG) = — (/?*(s) for every s€R. 

Hence, 

y?* (s + 2 In G) = -<p*(s + In G) = (f* (s) for every s ER, 

and the periodicity stated for (p* follows. According to the definitions in (16.91) and the 
expression of the right-hand interface yo, one obtains that 

V*(s + In GO = G- 5 e- 5 7o(yo - Ge s ) and <p*(s) = e~ 5s f (y - e s ). 
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Assume, for instance, that s is such that y — e s G (0, 1). Then, y — Ge s G (0, G) and, 
using the definition of fo in such an interval, we have that 

fo(y - Ge°) = -G 5 f { yj ^ - l) = -G 5 f (y ~ Ge*). 

The same argument applies in (0, G) and in all the intervals obtained in every step of this 
construction. This completes the proof. □ 

As we mentioned above, by reflection of the solutions constructed above, we obtain 
analogous results for solutions with the left-hand interface. We note that, in both cases 
(the left- and right-hand interfaces), the behavior of the self-similar profile F of (I2.5P close 
to the interfaces is described by means of /o(2/;G 2 ). The behavior corresponds to TW 
"travelling" in the opposite directions with A = +1 is given by the profile fo(y; G{). 

6.2. Uniqueness of oscillatory TW solutions. In Section m we have proved unique- 
ness of the solution to the Cauchy problem for the equation (12.21) with given initial con- 
ditions (14. 3p . when aj ^ at least for some values of j. It is clear, due to the regularity 
conditions at the interface y , that this analysis of uniqueness does not apply to the os- 
cillatory travelling wave solution at this point. In fact, it is not difficult to check that 
uniqueness fails and, for instance, the solutions / = 0, the positive solution constructed 
above, and its negative counterpart are three different solutions of the same problem with 
analogous conditions at the interface y = y^. 

In order to understand the structure of the set of solutions of this type, it would be 
interesting to fix an additional condition that provides uniqueness of the TW. We next 
prove that uniqueness holds by adding to the previous conditions at the interface yo, a 
fixed transversal zero ao of the solution and a sign to the function close to the fixed zero. 
Without loss of generality, we next assume ao = and positivity of the solution to the 
left-hand side of a . The result is stated as follows: 

Theorem 6.2. There exists a unique solution of the equation {\6.1\i with the interface yo 
with fv\yo) = for j = 0,1,2,3,4 and satisfying /(0) = and f positive in (— e, 0). 
Moreover, the solution is given by fo(y), with the algebraic construction above. 

Proof. Let f(y) be a solution of (16.11) satisfying the stated assumptions. Define the 
auxiliary function g(y) = f(y) + p(y), with p(y) a polynomial of the fourth degree, 

4 j 

p(y) = £ -Jr > 

i=l 

with rrij = (fo — /)^(0). By the definition of g and taking into account the properties of 
/ and f around a Q = 0, it is not difficult to check that g(y) satisfies the same equation 
as f on (s,e) and that g^(0) = /&')(o) for j = 0, 1, ...,4. 

Hence, it follows from the uniqueness result of the solutions to the initial value problem 
stated in Section HJ that g = fo at least in this interval, and, by a continuation argument 
(applied for instance to the integral version of the equation), we conclude the equivalence 
of both functions and its derivatives for every y < yo- 
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In particular, we have at the interface yo that g^\yo) = for j = 0, 1,2,3,4. Taking 
into account the definition of g and that fo and / satisfy also these regularity conditions 
at yo, it follows that p^\yo) = for j = 0, 1, 2, 3, 4, whence p(y) = and the uniqueness 
result follows. □ 

This allows us to establish a result concerning the structure of the set of oscillatory 
TW solutions to equation (16.11) with finite interfaces. 

Theorem 6.3. Assume that f(y) is an oscillatory TW solution of (\6.1\i with a finite 
interface and necessary regularity assumptions at the interface. Then, there exist real 
values 7 > and D such that 

f(y) = ± 1 5 fo( y - + D). 

Proof By using the rescaling properties of the solutions, it is enough to define an auxiliary 
function 

g(y) = ±\ 5 f(l + C) 

with A > and C such that y is the interface of g(y) and g(0) = 0. We also choose the 
sign in this definition in order to have all the hypotheses in Theorem 16.21 to hold, whence 
we deduce that g(y) = fo(y). It is not difficult to check that the result follows for 7 = j- 
and £> = -§■ □ 

Both results have strong implications concerning the properties and the character of the 
oscillations of the solutions and therefore, the behaviour of solutions close to the interface. 
In particular, we prove the following: 

Proposition 6.4. Assume that f(y) is an oscillatory TW solution with the interface yo. 
Then there exists a finite limit: 

lim sup ?~~)s = C. 

y^yo 

Proof. We prove the result for fo(y). Denote by y\ G ( — 1,0), the absolute maximum 
point of fo(y) in the interval (—1,0) and define 



y n +i = (y n + 1)G, with G 



yo 



1+2/0 ' 

From the properties of / an d its construction, it is not difficult to see that, so defined, 
y n are the absolute maximum point of \fo{y) \ in every interval I n of the extension defined 
in the algebraic construction of fo- Taking into account the definition of fo and y n yields, 
as n — > 00, 

_ G-(i- yi )G"+ 1 H 1/faOI _ G5("+i>/fa) r 
Vn ~ i=g allCl fe^F ~ d^-yn)5 °' 

so the result follows. □ 

As a straightforward consequence of the construction above, we also obtain the structure 
of the zeros of the solution f constructed: 

Corollary 6.5. Every zero of the solution fo{y), except the one corresponding to the 
interface, is transversal, i.e., f ^ 0. 
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7. Continuity geometric extension to m « 1 

As the next natural step, we consider equation (13.131) with a ~ 0. We then are going 
to perform the same geometric (but not algebraic!) construction explained in the steps 
(i)-(iii) in Section IBTTl where the only difference is that the basic profile fo(y) is not given 
explicitly as in (16.21) . We can use the standard continuous dependence of local solutions 
of the ODE (13.131) on the parameter a ~ 0. Here, we mean a class of good locally and 
uniquely extensible solutions without strong singularities and finite interface points. Of 
course, we crucially need the zero transversality property fixed in Corollary 16.51 Good 
local properties of such solutions are checked by standard contractivity techniques (as in 
Section l3T2"j) . so we arrive at a local continuity ("homotopy a-deformation" ) result: 

Theorem 7.1. Equation ( |5.^p with both signs " ± " (i.e., for A = =Fl) has a nontrivial 
periodic solution tp*(s) of changing sign for all a ~ 0. 

It is worth recalling here that, in a global sense, for a sufficiently far away from 0, peri- 
odic solutions do not exist; cf. Proposition 15. 21 and the heteroclinic bifurcation phenomena 
in (15331 . 

Notice that the existence result in Theorem 17.11 is not associated with the "hyper- 
bolicity" of the periodic orbit y?* in this range. Anyway, Theorem 17.11 altogether with 
Propositions 15.21 and 15.31 though not covering the whole a-range, provide us with rather 
solid mathematical evidence on existence of periodic solutions and also on the type of 
heteroclinic bifurcations, at which these disappear. 
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